Laser acceleration of ion bunches at the front surface of overdense plasmas 
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The acceleration of ions in the interaction of high intensity laser pulses with overdense plasmas 
is investigated with particle-in-cell simulations. For circular polarization of the laser pulses, high- 
density ion bunches moving into the plasma are generated at the laser-plasma interaction surface. 
A simple analytical model accounts for the numerical observations and provides scaling laws for the 
ion bunch energy and generation time as a function of pulse intensity and plasma density. 
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The study of the interactions between ultra-intense 
laser pulses and plasmas has proved to be a very rich soil 
where technological progress and fundamental physics 
meet each other. Particularly intriguing is the concept of 
laser-plasma based ion acceleration. From astrophysics 
0, to medical hadrontherapy , from proton radiog- 
raphy Q, to nuclear physics from proton imaging 
techniques Q, to nuclear fusion the problem of ac- 
celerating and manipulating charged particles with laser- 
plasma interactions offers a series of challenges ranging 
from fundamental to applied physics, thus a clear under- 
standing of the basic mechanisms is mandatory. Several 
recent experiments have reported the emission of ener- 
getic ions from solid targets 0. It is still a matter of 
debate whether the ions are mainly accelerated at the 
rear surface of the target (by the field generated by fast 
electrons escaping in vacuum ||) or at the front surface 
involving phenomena such as acceleration by a collision- 
less electrostatic shock 0, 0, , by a solitary wave 
or by ion trapping in a propagating double layer [T3 |. 

In this work we elucidate an even more basic process of 
ion acceleration in cold plasmas, purely related to the for- 
mation of an electrostatic field due to the action of the 
laser pondcromotive force on the electrons and, conse- 
quently, on the ions via space charge displacement. This 
investigation shows both the necessity of a kinetic de- 
scription of this process and the fundamental role played 
by the laser light polarization by showing the differences 
between circular and linear one. It will be shown by 
particle-in-cell (PIC) simulations that circularly polar- 
ized light gives rise to a "pulsed" acceleration and pro- 
duces ion bunches directed into the target. A simple 
analytical model is used to explain the acceleration dy- 
namics and for the deduction of scaling laws that relate 
the interaction parameters to the energy acquired by the 
ions. With respect to other concepts for laser ion accel- 
eration, the present mechanism with circularly polarized 
light leads to very high densities in the bunches, as might 
be of interest for problems of compression and accelera- 
tion of high-density matter. 

We consider a laser pulse impinging on a cold, step- 
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FIG. 1: Ion density ni (top, solid line) electrostatic field E x 
(top, dotted), and (x,p x ) phase space projection (bottom) 
at different times (labels) from ID PIC simulations. Run 
parameters are ai = 2, n e /n c = 5, I = 2Al- Eq = m e tuc/e. 



boundary, overdense plasma with no/n c = lo 2 /uj 2 l > 1, 
where uq is the initial electron density, n c = m e io\ /Aire 2 
is the critical density for a laser with carrier frequency 
lot,, <-<Jp is the plasma frequency and m e , e are the elec- 
tron mass and charge. The laser field amplitude will 
be given in units of the dimensionlcss parameter ul = 
(eEr,/m e ujr,c). In the PIC simulations, the resolution is 
high enough to resolve low density regions, sharp gra- 
dients and the dynamics of both electrons and ions by 
taking at least 20 particles per cell at n e = n c and a spa- 
tial resolution better than 0.1d p where d p = c/oj p . For 
reference, in all simulations shown the target boundary 
is located at x = and the laser impinges on it at t = 0. 

The essential features of the generation of ion bunches 
can be seen from one-dimensional (ID) simulations. The 
laser pulse is circularly polarized, incident from the left 
side and propagating along the x axis; its amplitude rises 
for 6Tl (where Tr, = 2-k/ujl is the laser period) up to the 
value ar, = 2 and then remains constant. The plasma 
has a slab profile with uq = 5n c and thickness t = 2\r,. 
Electrons and ions are assumed to be cold, T e = = 0, 
and ion charge and mass number are Z = 1, A = 1. 

Fig- U\ shows the ion density profiles and (x,p x ) phase 
space projections at different times. At t = 8Tl a sharp 
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FIG. 2: (Color online) a): Schematic of the model profiles 
of the ion density (green, thick) the electron density (blue, 
thin) and the electric field (red, dashed) at the initial stage 
when the electrons are in equilibrium and the ions have not 
moved yet. The label A indicates the electron front where 
the laser evanescence starts, B indicates the point where the 
ponderomotive force vanishes so that AB = l s . b): Compar- 
ison between the ion bunch velocity observed in simulations 
(data points) and the model prediction (lines), as a function 
of the laser pulse amplitude ah and the normalized density 
N = n i0 /n c . The cases studied are: N = 100, Z/A = 1/2 
(red, stars); N = 20, Z/A = 1 (blue, crosses); N = 5, 
Z/A =1/2 (green, diamonds). 



peak of the ion density reaching up to seven times the 
initial value is shown which moves inwards the plasma 
with an average velocity of about 0.013c leaving behind 
a low-density shelf with velocity linearly decreasing. At 
t = 12Tl, after the density at the peak has reached the 
maximum value of about 7no and the velocity has reached 
the value of 0.02c, a second density peak has appeared 
on the right of the first peak and moves into the plasma 
at the speed of about 0.04c with the (x,p x ) distribution 
now clearly a two-valued function of x. 

A simple model for this ion bunch generation can be 
inferred. When the laser pulse impinges on the plasma 
surface, electrons are quickly pushed inward by the pon- 
deromotive force, i.e. the steady part of the v x B force. 
(Notice that the oscillating part of the v x B force is zero 
for circular polarization; this is the reason why, as shown 
below, the interaction regime is completely different for 
linear polarization) . The electrons pile up leaving behind 
a charge depletion layer and giving rise to an electro- 
static field E x back-holding them. We assume that they 
quickly reach an equilibrium position where E x balances 
the ponderomotive force exactly. Let us take for sim- 
plicity a linear profile of E x both in the depletion layer 
(E x = Eox/d for < x < d) and in the compression 
region (E x = Eq[1 — (x — d)/l s ] for d < x < d + l s ), which 
implies a uniform electron density n P Q in this region, see 
Fig. El The parameters Eq, n p o, d and l s are related by 
the equations Eq = A-Ken^d (due to Poisson equation), 
no(d + l s ) = n p ol s (due to global charge conservation), 



and Eoenol s /2 ~ 2I L /c (due to the balance between the 
total radiation and electrostatic pressures). 

The electrostatic field on the ions starting at initial 
positions xq < d is a constant over the trajectories of the 
ions, E x = E x (xq), and increases with xq. Thus, these 
ions will never reach those with xo > d, the ion density 
will decrease and a ion "shelf" will be formed, the field 
on the leftmost ion layer being zero. As for the ions with 
initial position in the compression region d < xq < d+l s , 
assuming that the electrostatic field is a function of xo 
for these ions too is consistent with the assumption that 
the electrons follow a quasi-equilibrium dynamics and 
the ponderomotive force equals the electrostatic field at 
any time, the total pressure being always 21^/c. Being 
the electrostatic field a linear function of Xq , all ions will 
reach the point B at the same time and the ion density 
will assume an infinite value there, i.e. the hydrodynamic 
description breaks down. The ion bunch is thus due to 
accelerated ions that cross the point B and are injected 
into the unperturbed plasma region, where they can move 
ballistically provided that the charge unbalance is neu- 
tralized by electrons accompanying the bunch. The total 
number of accelerated ions per unit surface is ti p qI s . 

Note that the model assumes the point B, where the 
field vanishes, not to move during compression, and pre- 
dicts ions to have a a flat-top velocity distribution ex- 
tending between and v max and thus the density to 
be uniform in the region behind the ion front. Ac- 
tually, since for the evanescence length l s we expect 

— 1/2 

l s ~ d p ~ n e , during the compression of the ion fluid 
the field will tend to penetrate deeper into the plasma 
keeping the field at the surface and the total electrostatic 
pressure constant. Thus, ions beyond point A will be ac- 
celerated by a field decreasing in time and will get to the 
breaking point later. This effect causes the ion bunch to 
be more localized both in coordinate and velocity space. 

Our model thus gives a scenario in qualitative agree- 
ment with the numerical observations and also provides 
quantitative estimates and scaling laws. Denoting as 
= Co + d the initial position of an ion with mass mi 
and charge in the compression region, (0 < Co < h), 
the force acting on this ion is given by Fi = qiEo(l — 
Co /l s )- Thus, The velocity of an initially immobile ion 
is Vi = (qiEo/n%i){l — Co/h)t and the position is Xj = 
(qiEo/2m,i)(l — (o/l s )t 2 + xq. The "breaking" time at 
which all ions reach point B is r, = \j2l s m,ij 'qiE^. The 
maximum velocity is thus v max = y / 2l s qiEo/mi = 2v a , 
being v a = I s /t% the average ion velocity, i.e. the velocity 
of the laser front. By relating the model parameters to 
the laser intensity and plasma density we obtain 
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where l s ~ d p has been assumed. The predicted velocity 
scaling has been tested via numerical simulations, yield- 
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FIG. 3: a) Ion bunch acceleration from a thin {I — Q.IXl), 
solid-density (no/n c — 100) target irradiated by a short laser 
pulse with amplitude a = 5 and duration of QTl (FWHM). 
The ion density profile (thick line) and the (x,p x ) phase 
space distribution (dotted) are shown at t — 13Tl- b) Ion 
density (top) and phase space (x,p x ) projection (bottom) at 
t = 86Tl , for the same simulation parameters of Fig. but 
linear polarization of the laser pulse. 



ing a rather good agreement with the model, see Fig. 
12 For instance, for the parameters of Fig. the model 
predicts v a /c = 0.021 and Ti ~ 3.4T^; in the simulation 
the breaking of the ion front occurs at t k, 102l which, 
taking into account the pulse rise time of 6Tl, is in fair 
agreement with the model prediction. 

At the "breaking" of the ion profile with the formation 
of the bunch, the equilibrium between the electrostatic 
and the ponderomotive forces on electrons is lost and, if 
the laser pulse is not over, electrons rearrange themselves 
to provide a new equilibrium. The process of bunch for- 
mation might then restart although with different initial 
conditions, so that the simple model is not adequate any- 
more for quantitative predictions. In Fig.^a second wide 
bunch actually appears at t ~ 15Tl and then propagates 
into the plasma. At later times, bunches with shorter 
duration are generated at a higher rate. 

For potential applications, it is interesting to study the 
possibility to produce a single bunch from a thin solid tar- 
get. Results for the case of an n /n c — 100, Z/A =1/2 
target with a thickness I = 0.1 Xl and a short pulse with 
intensity cll = 5 are shown in Fig. [3] The temporal en- 
velope for the field is given by sin 2 (7rt/r) with r = 12. A 
single bunch of ions with velocity w 0.017c is produced at 
t « 5.57l and leaves the target from the back side. The 
ion bunch charge is completely neutralized by electrons 
both inside and outside the target. 

Collisional ion stopping (not included in the PIC code) 
should not affect ion bunch generation strongly for thin, 
low-Z targets. In fact, taking for the example of Fig[3| 
A = 0.5 /im, so that no = 4 x 10 23 cm~ 3 , the stop ping 
power in the cold plasma given by the Bohr formula [l4| 
is ~ 0A9Z 2 MeV//im for ions with v = 0.017c and charge 
Z = A/2, and the stopping length is thus w lAZ^ 1 fixa, 
lower than the target thickness in the example if Z < 28. 
For lower densities, typical of preformed plasmas or, e.g., 
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FIG. 4: (Color online) Contours of ion density rii/n c from 2D 
simulations at three different times. The laser pulse impinges 
from the left and has a transverse FWHM of 4Al- Notice the 
different scales on the x and y axis. 



foam targets, collisional ion stopping is even less impor- 
tant. Collective ion stopping in the plasma is included in 
the PIC simulation, and is found not to affect ion trans- 
port in the regimes addressed by our simulations. 

The mechanism of ion bunch formation outlined above 
works cleanly for circular laser polarization. For lin- 
ear polarization the oscillating part of the v x B force 
causes strong electron heating (up to MeV energies) 
which changes the regime of ion acceleration qualita- 
tively. Fig. |3 shows results from a simulation with pa- 
rameters identical to Fig.^but the laser polarization that 
was changed to linear. No clear propagating ion bunch as 
that of Fig. is observed. The hot electrons form a De- 
bye sheath leading to ion acceleration towards vacuum at 
the back surface as soon as they reach it while a few ions 
are also accelerated from the front surface towards vac- 
uum in a similar way (Fig.|3J. Near the laser pulse front, 
the phase space structures look similar to those gener- 
ated by the shock acceleration investigated in Ref . [lOj at 
higher intensity. For the present parameters, sheath ac- 
celeration dominates at the time shown. In contrast, for 
circular polarization, it is found that most of the elec- 
trons are heated up to energies of a few keV behind the 
ion bunch and at the laser-plasma interaction surface. 

Notice that in the case of shock acceleration (see e.g. 
Ref. [HI) the laser pulse acting as a piston drives a shock 
wave into the plasma which reflects the ions thus acceler- 
ated up to velocities twice the shock speed. In the present 
case the fastest ions have twice the piston velocity, i.e. 
the velocity of the laser reflection front at the breaking 
time, and they come from behind the front. It is also 
worth to notice that the present mechanism of ion bunch 
formation is of electrostatic and kinetic nature while a 
purely hydrodynamic description is not adequate. 
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FIG. 5: (Color online) The (y,Px) and (p x ,p y ) phase space 
projections [a)-b)] and the energy spectrum [c)] of ions from 
the 2D simulation of Fig.gl at t = 47T L . 

Wc performed 2D simulations to show that the mecha- 
nism for bunch formation is also at play in more than one 
dimension, and in particular in the case of a laser pulse 
with a finite transverse width that bends the plasma sur- 
face (hole boring). Fig0] shows results for a case where 
the laser pulse has a Gaussian intensity profile in the 
transverse (y) direction with halfwidth w = 4Al, dura- 
tion t = 12Tl, and peak intensity ar, = 2. The target 
thickness is I ~ 0.5Al and the density is = hn c . At 
t = 18Tl, we see the gaussian-shape compression front 
produced by hole boring. At t — 301l, the ion bunch 
has formed and propagates into the plasma, keeping the 
shape of the compression front. At t = 427Y, the bunch 
ions in the central region have left the target and prop- 
agate in vacuum. Another run was performed with the 
same parameters but a "tighter" pulse, i.e. w = 2A^. 
and almost identical results were found. 

FigEl shows the ion phase space projections and the 
energy spectrum for the simulation of FigQ] at t = 47Tl. 
The dependence of the longitudinal momentum p x on y 
[FigEl a)] and the energy spectrum [FigGO c)] resemble 
the intensity modulation of the laser pulse. The relative 
spread in the transverse momentum [FigEl b)] is very 
small, (ip 2 y ) I (pl)) 1/2 ~2x 1CT 2 , showing that also in 2D 
the acceleration is perpendicular to the target surface. 
In the w = 2Xl case we find slightly slower ion velocities 
and almost the same beam divergence. 

In some runs, some "rippling" of the ion density front is 
observed. This rippling appears to be due to electron os- 
cillations at the surface (around their quasi equilibrium 
position) excited due to the finite rise time and the curva- 



ture of the laser pulse. The rippling is "imprinted" in the 
ion density at later times, but has no noticeable effect on 
the dynamics of ion acceleration. We conclude that 2D 
effects do not qualitatively affect ion bunch formation. 

In conclusion, a regime of laser ion acceleration in cold 
overdense plasmas by circularly polarized pulses has been 
characterized by means of PIC simulations and a simple 
model. The dynamics is qualitatively different from the 
case of linear polarization where ion acceleration is driven 
by fast electrons. With the present scheme one may ob- 
tain moderate ion energies but very high ion densities 
and low beam divergence. 
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